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For Phys.Rev.B 

It is shown here that Kohn-Sham equations cannot be 
derived from Hohenberg-Kohn theory without an additional 
postulate. Assuming that a functional derivative with respect 
to total electron density exists leads in general to a theory 
inconsistent with the exclusion principle. A mathematically 
and physically correct variational theory of the Kohn-Sham 
model can be developed using functional derivatives with re- 
spect to orbital densities. These partial (Gateaux) deriva- 
tives can be constructed explicitly from general N-electron 
theory and are defined throughout the orbital Hilbert space. 
This theory is consistent with the local density approxima- 
tion (LDA), but does not in general imply multiplicative local 
exchange-correlation potentials. Progress beyond the LDA in 
condensed-matter physics requires development of methodol- 
ogy for nonlocal exchange and correlation potentials. 

71.10.-|-x,31.10.-f-z 



I. INTRODUCTION 

Hohenberg and Kohn^ proved, for nondegenerate 
states, that the ground-state total electronic density 
function determines the external potential function act- 
ing on an interacting N-electron system. This is easily 
extended to a spin-dependent density, so spin indices 
and summations are assumed here, but omitted from 
the notation. Implementation of the implied density- 
functional theory (DFT) followed only after introduction 
of the orbital model of Kohn and Sham^, who applied 
the Hohenberg-Kohn theory of the universal ground-state 
functional Fg = E — V to the kinetic energy functional 
Tg, to which Fg reduces for nonintcracting electrons. For 
noninteracting electrons, the density function p(r) is a 
sum of orbital densities pi = (t>*4>i, one for each occupied 
orbital function of the Kohn-Sham model state. DFT 
is developed in many review articles and monographs, 
exemplified by^'*. 

It is widely assumed that Hohenberg-Kohn theory im- 
plies the existence of an exact theory that can be rep- 
resented by the Kohn-Sham model, using only multi- 
plicative local potential functions. It is shown here 
that this assumption is based on a misconception of 
the implications of Hohenberg-Kohn theory. If the ac- 
cepted interpretation were correct, an exact Thomas- 
Fermi theory (TFT) would exist, in conflict with Kohn- 
Sham DFT and with the exclusion principle for elec- 
trons. The assumed existence of exact local exchange- 



correlation potentials leads to several problems or para- 
doxes, resolved by the present analysis. Recent quanti- 
tative tests^ confirm a longstanding result of optimized 
effective potential (OFF) calculations, that the best pos- 
sible local exchange potential does not reproduce unre- 
stricted Hartree-Fock (UHF) ground-state energies and 
densities^. The exchange-only limit of DFT linear re- 
sponse theory'' is inconsistent with the time-dependent 
Hartree-Fock theory of Dirac®, due to failure of locality^. 

The present analysis has important practical impli- 
cations for computational methodology relevant to the 
electronic structure of large molecules and solids. In a 
formally exact theory, deduced directly from N-electron 
theory, nonlocal exchange-correlation potentials can be 
constructed that are free of electronic self-interaction and 
not inherently restricted to weak electronic correlation 
or to short-range interactions^". Progress beyond the 
local density approximation requires the timely develop- 
ment of computational methods appropriate to nonlocal 
potentials'^. 



II. DFT PRACTICE VERSUS DFT CONCEPTS 

Practical applications of DFT use the Kohn-Sham 
Ansatz, p — riipi, as a postulate for the density. Oc- 
cupation numbers are unity for occupied orbitals. Ki- 
netic energy is defined (in the nonrelativistic limit) by 
the Schrodinger functional of occupied orbital functions, 
T ~ where t = — ^V^ in Hartree units. 

Given an external potential function w(r), the potential 
energy V = X^i'^iC*!^!*) orbital functional, as is 

the Hartree energy E^. By the Kohn-Sham Ansatz, any 
assumed Exc[p] is a functional of the occupied orbital 
functions. Given F = T + Eh + E^c as an orbital func- 
tional, derivation of the orbital Euler-Lagrange (OEL) 
equations 



= {e, - v{r)}(j)i 



(1) 



is elementary, using Lagrange multipliers for indepen- 
dent orbital normalization — l,i < N . Linear oper- 
ator J- is defined by the orbital functional derivative 



SF 



(2) 



F is assumed to be defined such that J- is hermitian. 
For noninteracting electrons, and also in the local-density 
approximation (LDA) , when E^^ is the integral of a local 
function of p, the OEL equations (1) reduce explicitly to 
Kohn-Sham equations^ with local potentials. 
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Eqs.(l) imply equivalent equations^^ 



5Fs = j d^rvF{r)Sp{r), 



(4) 



VFi{r) 



v{r), 



(3) 



to be solved for the orbital density functions p;. The 
Lagrange multipliers Cj are to be determined so that in- 
dependent normalization conditions = J Picfr = 1 
are satisfied. These modified Thomas-Fermi equations 
are operationally equivalent to the Kohn-Sham or OEL 
equations, consistent with the exclusion principle. 

This practical derivation makes no use of Hohenberg- 
Kohn theory except to justify particular approximations 
to E^c- The density function is constructed after the fact 
from solutions of the variational equations. To define the 
model state, some rule ^ — > $ must be postulated to 
supplement the Kohn-Sham Ansatz for the density. Ex- 
amples are p<j> = (Kohn-Sham), \1/ = $ (unrestricted 
Hartree-Fock), and (^^l^*) = 0, the Brueckner-Brenig 
condition that selects the model state of greatest weight 
in $1'^'^^. A derivation starting from Hohenberg-Kohn 
theory is more problematic, since a formal theory of al- 
lowable densities is required^''. In the orbital theory de- 
scribed above, this is avoided by recognizing that densi- 
ties freely constructed from orbital functions in the usual 
Hilbert space suffice for a valid N-electron theory. 



III. INCOMPLETENESS OF 
HOHENBERG-KOHN THEORY 

The Hohenberg-Kohn ground-state functional Fg is de- 
fined only for normalized total densities, J p = N. Stan- 
dard variational theory constrains normalization using 
the modified Lagrange functional Fs + V — p[J p — N], 
which is to made stationary by varying p with uncon- 
strained normalization^^, p is adjusted to satisfy the 
normalization constraint. This requires the definition of 
Fs to be extended into all infinitesimal neighborhoods of 
normalized densities consistent with the orbital Hilbert 
space. Without such an extended definition, Hohenberg- 
Kohn theory is incomplete: it cannot determine the vari- 
ational Euler-Lagrange equations for p. 

This point is evident in the rigorous derivation of En- 
glisch and Engiisch^^. In their Eq.(4.1) they propose 
to complete their derivation by inserting a single con- 
stant, which would play the role of the usual chemical 
potential p if the implied Euler-Lagrange equations were 
valid. However, comparison with Eqs.(3), implied by 
physically correct Kohn-Sham equations, indicates that 
a single constant does not suffice to enforce the exclusion 
principle^*. Independent normalization of each orbital 
density is required. The derivation iu"'^^, restricted to nor- 
malized total densities, cannot distinguish between one or 
several constants, all of which drop out of the variational 
integrals. As shown here, this is a crucial limitation of 
Hohenberg-Kohn theory. 

If the functional differential 



is uniquely defined for unrestricted variations, it deter- 
mines a total (Frechet^^) density functional derivative 



SFs 
6p{r) 



VF{r). 



(5) 



This is stated as the unique definition of a density func- 
tional derivative in standard DFT literature'^'^. If such a 
Frechet derivative exists, varying the Lagrange functional 
using standard variational theory-^^ implies the Euler- 
Lagrange equation 



vpir) = p-v{r). 



(6) 



This is a Thomas-Fermi equation, to be solved by vary- 
ing a trial density for specified p, then adjusting p to sat- 
isfy the normalization condition. It is inconsistent with 
Eqs.(3) unless all are equal^*'^". 

Because Thomas- Fermi theory (TFT) enforces only to- 
tal density normalization, it can determine only a single 
Lagrange multiplier p. This does not exclude solutions 
that violate the exclusion principle. The example of the 
lowest ls2s state of an atom containing two noninter- 
acting electrons is discussed below. In this example, two 
independent orbital energies eis and €2s are required, not 
just the single parameter p. The conclusion is that the 
required Frechet density functional derivative cannot ex- 
ist in a physically correct theory for any compact system 
with more than one electron of each spin. 



IV. THE EXCLUSION PRINCIPLE IN DFT 

In any independent-electron orbital model, the exclu- 
sion principle requires independent normalization of each 
orbital density. This is explicit in any physically correct 
theory of noninteracting electrons, implicit in the Kohn- 
Sham equations, and is a standard aspect of Hartree-Fock 
equations. If the Kohn-Sham density decomposition is 
added as a postulate to the analysis of Englisch and En- 
glisch, the constraint of independent orbital normaliza- 
tion requires independent variation of the orbital partial 
densities. Such density variations, driven by free (unnor- 
malized) variations of occupied orbital functions, deter- 
mine Eqs.(3) as Euler-Lagrange equations for the orbital 
densities^^. These equations are expressed in terms of 
partial (Gateaux^^) functional derivatives of the univer- 
sal functional F=E-V. Independent normalization of each 
orbital density introduces independent Lagrange multi- 
pliers, which are the one-electron energies in the corre- 
sponding orbital Eqs.(l). 

Eqs.(3) and (6) can be reconciled only if all orbital 
energy eigenvalues are equal^^'^°. This is seen here to 
be a direct consequence of the different normalization 
constraints. That independent normalization of the or- 
bital densities is physically correct, specifically enforcing 
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the exclusion principle in Hartrce-Fock theory and Kohn- 
Sham DFT, is obvious in the example of an atom with 
two nonintcracting electrons of parallel spin. The ground 
state is ls2s ^S. The density constraints are: 

DFT : j pu = l, J P2s = 1 

TFT: j{pu+P2s)=2. 

This TFT constraint does not exclude the nonphysical 
solution 

TFT : J pu=2, J p2s = 0, 
which violates the exclusion principle. 

V. GATEAUX FUNCTIONAL DERIVATIVES 

Euler-Lagrangc equations arc implied for interacting 
and nonintcracting electrons, respectively, only if the 
definitions of F^ and Ts can be extended to unnormal- 
izcd densities generated by infinitesimal variations un- 
constrained in the orbital Hilbert space. Given F^ as a 
functional of total density, it is also an orbital functional 
as a result of the Kohn-Sham Ansatz. The orbital func- 
tional derivative in the OEL Eqs.(l) exists if 



(7) 



defines a unique functional of the occupied orbital func- 
tions. For variations about a stationary state, Eqs.(l) 
imply that 



SFs = J2'^i / (fr{e^-v{r)}6p,{r), 



(8) 



since Spi = {S(j)*(j)i + cc}. This proves, by construction, 
the existence of a Gateaux functional derivative whose 
value for such variations is 



SF, 
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■v{r), 



(9) 



for i < N. A total (Frechet) derivative is defined only 
if these partial functional derivatives are all equal, which 
requires all to be equal. 

Gateaux functional derivatives exist in general 
throughout the orbital Hilbert space. They can be con- 
structed from any well-defined density functional, using 
the Kohn-Sham orbital decomposition of the total model 
density function. They define local potentials indexed by 
orbital subshells, in general not a global local potential 
unless the orbital index drops out. For example, for any 
orbital functional the functional differential 



SF = J2ni f d^v{5(t>*J^(j)i + cc} 

i 



is uniquely defined tov i < N throughout the orbital 
Hilbert space if the orbital functional derivative ^c^, ex- 
ists. This determines the Gateaux derivative 
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(t>*(t>, 



= VFi{r), 



(11) 



valid for arbitrary variations in the orbital Hilbert space. 
This logic validates Eqs.(3) in a density-based theory. 
Because this orbital-indexed multiplicative local poten- 
tial is singular at any node of (f)i, these Gateaux deriva- 
tives can only be equal for all index values if the nodes 
exactly cancel. A Frechet derivative can occur only in 
such a special case, which is very difficult to construct 
unless the functional is defined as in the local density 
approximation. 

Hence extending the definition of Fs[p] so as to de- 
termine TFT cqiiations conflicts in general with the ex- 
clusion principle. In contrast, an extended functional 
Fs[{pi}] and its Gateaux derivatives exist, implying 
Eqs.(3). The corresponding OEL equations reduce to 
Kohn-Sham equations in the LDA, but locality of ki- 
netic, exchange, and correlation potentials is not implied 
for more general energy functionals. The mathemati- 
cal statement of this situation is that the Hohenberg- 
Kohn argument does not imply the existence of a Frechet 
functional derivative, equivalent to a multiplicative lo- 
cal potential function, for any system with more than 
one Fcrmi-Dirac electron of each spin. This conclusion 
does not conflict with the rigorous theory of Englisch 
and Englisch^^. Their argument, restricted to normalized 
ground-state total densities, cannot distinguish between 
Frechet and Gateaux functional derivatives. 



VI. ORBITAL FUNCTIONALS IN AN EXACT 
THEORY 

Formally exact orbital functionals for kinetic, ex- 
change, and correlation energy can be deduced from stan- 
dard N-electron theory. They can be incorporated into 
a formally exact mean-field theory in which the OEL 
equations extend Kohn-Sham theory to nonlocal poten- 
tials (orbital- indexed local potentials). Progress beyond 
the local-density approximation requires development of 
computational methodology for such nonlocal potentials. 

The variational theory of independent-electron models 
is most simply developed as an orbital functional the- 
ory (OFT)^°. For an N-electron eigenstate such that 
{H — E)^ = and any rule ^ ^ that determines a 
model or reference state unsymmetric normalization 
($1*) = ($1$) = 1 implies that E = ($|if|*) = Eq+Ec, 
where Eq = {^\H\^) is an explicit orbital functional. 
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and Ec = — defines the correlation energy. 

Restricting the discussion to nondegenerate states, $ 
is a Slater determinant constructed from occupied or- 
bital functions {4>i}, with occupation numbers rij = 1. 
Spin indices and sums are assumed here but are sup- 
pressed in the notation. For orthonormal orbital func- 
tions, Eq = T + U + V, where 



T = ^ni{i\i\iy, V = J2Mi\v\i), 



(12) 



= VFi{r) =ei-v{r), i = l,---,N. 



(19) 



These equations reproduce Eqs.(3), derived directly from 
the OEL equations. This derivation generalizes Eq.(6) of 
TFT to incorporate the exclusion principlc^^. Orthogo- 
nality constraints consistent with independent variation 
of orbital densities follow from the operationally equiva- 
lent OEL equations. 



for t = — iy^. Two-electron functionals are defined for 
u = l/ri2 hy U = Eh + E^, where 

Eh = ^ ^n»%(ii|uKi); E^ = ~^ninj{ij\u\ji). (13) 

If Q = / - $$t^ = ($|if|* - $) = ($|iJ|e*). This 
implies an exact but implicit orbital functional^'' 



Ec = -{^H[Q{H -Eo~Ec~ iv)Qr'H\<S>) 



(14) 



for 7] — » 0-I-. Because Ec is determined by the antisym- 
metric model function $, its functional derivatives cannot 
introduce spurious electronic self-interaction. 

Orbital Euler-Lagrange equations follow immediately 
from standard variational theory^^, in terms of the or- 
bital functional derivatives 



5V 

niS(f>* 



5T . 5U 

tcpi] ^ = u(pi; 



v{v)(l)i; 



SEc' 
ni5(f>* 



Vc(i>i, 



(15) 



using u = Vh{r) + Vx, where Vh is the classical Coulomb 
potential, and Vx is the Fock exchange operator. Defining 
F = E — V = T + U + Ec, its orbital functional derivative 
is 
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T^i = {i + u + Vc}4>i. 



(16) 



For independent subshell normalization, Lagrange terms 
ei[(i|i) — 1] are subtracted from the energy functional. 
The variational condition is 



j ST^ni{5(l)*{T + v-ei}(l}i + cc) = Q, 



(17) 
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for unrestricted orbital variations in the usual Hilbert 
space. This implies the OEL equations 



= {ei - v}(t>i, i = l,---,N. 



(18) 



This theory can be expressed in terms of orbital den- 
sities, using the modified Lagrange functional F -\-V — 
n,ei[/ Pi — 1]. Using Eq.(ll), the stationary condi- 
tion is expressed as a relationship between local potential 
functions, 
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